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We study the three-nucleon force in chiral effective field theory with explicit A-resonance degrees 
of freedom. We show that up to next-to-next-to-leading order, the only contribution to the isospin 
symmetric three-nucleon force involving the spin-3/2 degrees of freedom is given by the two-pion- 
exchange diagram with an intermediate delta, frequently called the Fujita-Miyazawa force. We 
also analyze the leading isospin-breaking corrections due to the delta. For that, we give the first 
quantitative analysis of the delta quartet mass splittings in chiral effective field theory including the 
leading electromagnetic corrections. The charge-symmetry breaking three-nucleon force due to an 
intermediate delta excitation is small, of the order of a few keV. 

PACS numbers: 13.75.Cs,21.30.-x 



I. INTRODUCTION AND SUMMARY 



The importance of the A-excitation in the three-nucleon force (3NF) has been already realized fifty years ago by 
Fujita and Miyazawa in their seminal work [![. Their paper has been the seed for many meson-theoretical approaches 
i " i ' to the 3NF like the families of Tucson-Melbourne 0, Q , Brazilian Q or Urbana-Illinois 0, 0] 3NFs, see the review 
article Q and the recent general introduction Q. Nowadays, the appropriate tool to analyze the forces between 
' nucleons is chiral effective field theory, a program started by Weinberg |9|. There has been quite a sizeable body 
of further work on the structure of 3NFs within the framework of EFT, with appropriate references given below. 
Still, in the EFT with explicit deltas most investigations so far have considered the effects based on the leading 
pion-nucleon-delta coupling. In this work, we want to go one step further. In the two-nucleon system, we had already 
considered the effects of subleading ttNA couplings on the description of the peripheral phase shifts [1JJ . Such effects 
appear at next-to-next-to-leading order (NNLO) in the Weinberg counting. It is therefore appropriate to extend these 
considerations to the 3NF, more precisely to the A-contributions to the isospin-symmetric 3NF at NNLO and the 
leading isospin-breaking corrections to the 3NF. An important ingredient in latter type of forces stems from the mass 
splittings in the quartet of A-states (A ++ , A + , A , A"). These splittings receive contributions from the strong and 
■ electromagnetic interactions, thus one needs to consider an appropriate extension of the power counting including the 
. |_J ' explicit soft and hard photon effects. 

x : 

• The pertinent results of this investigation can be summarized as follows: 

a ■ 

_ _ i 

i) We have evaluated the A-contributions to the three-nucleon force up to NNLO and shown that the only non- 
vanishing topology is the two-pion exchange diagram with an intermediate delta resonance, commonly called the 
Fujita-Miyazawa force. This implies that the leading contributions to the l7r-4iV-contact and the 6iV-nucleon 
contact topologies are not saturated by the A. 

ii) We have compared the 27r-exchange 3NF in the EFT with and without explicit A degrees of freedom. We show 
that these representations lead to comparable results for the strength of the 3NF if the same data basis for 
pion-nucleon scattering is used to pin down the LECs in the two versions of the theory, see table |TJ 
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iii) We have analyzed for the first time the mass splitting in the delta quartet from strong and electromagnetic 
contributions with EFT up-to-and-including second chiral order. Note that chiral extrapolation formulae for 
the strong splitting have already been given earlier in Ref. [13]. The delta mass splittings are parameterized 
in terms of two independent parameters, cf. Eq. (|3.7[) . These can only be determined with large uncertainties 
since the available information on the various delta masses is fairly scarce and uncertain, see tables HTlIIII 

iv) We have shown how the proton-neutron mass difference can be eliminated from the delta-full EFT by field 
redefinitions, extending the method developed in (25[. This facilitates the calculation of the isospin-breaking 
effects to the 3NF considerably 

v) We have worked out the leading isospin-breaking contributions to the 3NF due to an intermediate A-excitation, 
see Fig. [1] The leading charge-symmetry conserving and charge-symmetry breaking contributions to the 3NF 
are given in Eq. (|3.24[) in momentum space. We also give the coordinate space representation. We estimate the 
contribution from the charge-symmetry breaking force to the 3N binding energy to be of the order of a few keV. 

The manuscript is organized as follows. In sec. [IT] we investigate the delta contributions to the 3NF up-to-and- 
including next-to-next-to-leading order in the so-called small scale expansion (SSE) [l(j. The calculation of the 
leading isospin-violating contributions is presented in sec. IIII1 As a first step, we analyze the mass splittings within 
the delta quartet and calculate the strong and electromagnetic contributions to the various particles, see sec. IIII Al 
In sec. IIII Bl wc discuss field redefinitions to eliminate the proton-to-neutron mass shift from the effective Lagrangian 
which considerably simplifies the calculation of the isospin-breaking effects. The isospin-breaking 3NFs due to explicit 
deltas are then worked out in sec. IIII CI 



II. A-CONTRIBUTIONS TO THE THREE-NUCLEON FORCE UP TO NNLO 

Our calculations are based on Weinberg's power counting Q utilizing the small scale expansion [To| . In this framework, 
irreducible diagrams with two or more nucleons which give rise to the nuclear forces are ordered according to the 
power v of the expansion parameter Q/A, where Q collectively denotes small pion four-momenta, the pion mass, 
baryon three-momenta and the nucleon-delta mass splitting and A is the pertinent hard scale. For an irreducible 
iV-nucleon diagram, the power v is given by: 

v = -2 + 2N + 2(L-C) + ^2ViAi. (2.1) 

i 

Here, L, C and Vi refer to the number of loops, separately connected pieces and vertices of type i, respectively. 
Further, the vertex dimension Aj is given by 

A, = dk + h,-2, (2.2) 

where hi is the number of baryon field operators and di is the number of derivative, insertions of M T and/or the 
delta-nucleon mass splitting, A = to a — tun- 1 

For the calculation of the isospin-conserving three-nucleon force (3NF) up to NNLO we use the effective chiral 
Lagrangian already given in [ll| . The only additional terms that need to be considered are the leading-order N N — > 
NA and NA — > JVJV contact interactions 

TfNNS^N + h. c. , (2.3) 

where N denotes the large component of the nucleon field, is the large component of the delta field, with i an 
isospin and /i a Lorentz index. Furthermore, t 1 and are Pauli isospin matrices and denotes the covariant spin 
vector. For more details on the notation and the effective Lagrangian the reader is referred to [H}. The contact 



1 Notice that we use the symbol A for both the spin-3/2 field and the NA mass splitting. It is, however, always evident from the context 
what is meant. 



(a) (b) (c) (d) 

FIG. 1: Leading isospin-conserving (diagram (a)) and isospin-breaking (diagrams (b-d)) contributions to the 3NF with a 
A-excitation. Solid, dashed and double lines represent nucleons, pions and deltas, respectively. Solid dots denote the leading 
isospin-invariant vertices, crossed circles are the isospin-breaking vertices from the S , (7(2)v-rotated effective Lagrangian as 
explained in section IIII Al 



interactions in Eq. (|2.3[) were already considered by van Kolck [l4j who worked out the corresponding contributions to 
the 3NF. As pointed out in Ref. |15j, matrix elements of the resulting 3NFs between antisymmetrized \NNN) states 
vanish, see also [l6[ for a related discussion. The reason for vanishing of these 3NF contributions can be understood 
already at the level of the effective Lagrangian. To that aim, let us rewrite Eq. (|2.3[) including the spin and isospin 

3/2 projectors P^ 2 and t^- 2 explicitly: 



T?NNS ll T i N = T"& 2 P3 /2 NNS^t 1 N, 



(2.4) 



with 



£ 3/2 - -S- 
3 — 3 



p3/2 = 



9^u ~ VpVv 



1 



(2.5) 



Here, denotes the baryon four-velocity and d the number of space-time dimensions. In four dimensions with 
V = (1, 0, 0, 0), the spin 3/2 projector reduces to 



P, 



3/2 



3*0 



(2.6) 



with a k denoting the Pauli spin matrices. It is straightforward to see that the antisymmetrized Feynman rules for 
the particular contact interactions vanish: 



d f2 ^)P 3 J 2 (1)44^2 = d /2 (2)P 3 kl /2 (2)r(a[A 12 = 0, 



(2.7) 



where the subscripts of the Pauli spin and isospin matrices refer to the nucleon labels and the antisymmetrization 
operator for momentum-independent interactions A12 has the form 



A12 = 



1-P 



12 



with Pi 2 = 



1 + CTi • a 2 1 + T\ • r 2 



(2.8) 



3/2 3/2 

The labeled projector operators £- ■ (X) and Pj (X) in Eq. (|2.7[) are defined according to 



lj ^ ' ~ 3 j ~ 3 j x ' 



p-/ 2 (x) = 



2 X _ * k 



xe{i,2}. 



(2.9) 



We are now in the position to discuss the 3NF contributions due to intermediate A excitations. Due to vanishing 
of the lowest-order contact interaction, the complete effect due the A is given by a single two-pion (27t) exchange 
diagram (a) in Fig. [TJ It is well known that the corresponding 3NF contribution is exactly reproduced at NNLO in 
EFT without explicit A-fields via resonance saturation of the low-energy constants (LECs) 03^4 accompanying the 
subleading mrNN vertices [17J: C3 = — 2c4 = — Ah\/(9A). Here, Ha denotes the leading ANir axial-vector coupling 
constant. 

Let us now regard the NNLO contributions. First, we emphasize that there are still no diagrams involving contact 
interactions at this order. In fact, since NNNA contact interactions contain at least two derivatives, they start to 
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Q 3 , no A 


Q 2 with A, fit 1 


Q 2 with A, fit 2 


Q 3 with A, fits 1,2 


a 


-0.70 








-0.70 


b 


-2.34 


-1.70 


-1.03 


-2.18 


d 


-0.89 


-0.42 


-0.26 


-0.83 



TABLE I: Values of the coefficients a, b and d entering the 2-7r-exchange 3NF in Eq. (|2.11[) in units of 1 for a and M w 3 for 
b and d. 



contribute to the 3NF only at N 3 LO. Therefore, the only diagrams which need to be considered are the ones which 
result from graph (a) in Fig. [1] by substituting one of the leading ttN A vertices by the subleading one which is 
proportional to the combination of LECs (63 + &§), see Refs. [TJ UM- 2 Since this vertex involves a time derivative, 
the corresponding 3NF contribution is shifted to higher orders due to the 1/mjv-suppression. We, therefore, end up 
with the conclusion that the 27r-exchange 3NF from diagram (a) in Fig. [T] represents the only additional contribution 
which arises in EFT with explicit A up to NNLO. 

Before ending this section, it is interesting to compare the strength of various pieces in the 27r-exchange 3NF in EFT 
with and without explicit A degrees of freedom which has the form 



where ql = pi ' — pi\ pi (pi ') is the initial (final) momentum of the nucleon i and 

Ff = 6 a ? [-a + bq l - q 3 - C (qf + qf)] - d e a ^ $ X % ■ a k . (2.11) 

The coefficients a, b and d can be expressed in terms of various LECs while c = 0, see Ref. [l8| for more details. In 
EFT without explicit A, the chiral expansion for these coefficients starts at NNLO (i.e. at order Q 3 ) where one has: 

7T 7T IX 

Here the superscripts refer to the chiral order. In EFT with explicit A fields, the dominant contributions arise already 
at NLO 

a (2) _ n l(2) %JA_ d (2) 2^_ 

with the corrections at order Q 3 being still given by Eq. (|2.12p . In Table HI we compare the values for the coefficients 
a, b and d in the theory with and without explicit A based on our determination [11] of the LECs Cj and 63 + b s . We 
remind the reader that fits 1 and 2 are based on the different values for the ttNA coupling constant \ia used as an 
input: the SU(4)/large-N c value h A = ig A /{2y/2) = 1.34 in fit 1 versus h A = 1.05 [H in fit 2. Clearly, at NLO this 
uncertainty in the value of h A directly transforms into the uncertainty in the coefficients b and d, see the third and 
fourth columns in Tabled At NNLO, however, both fits 1 and 2 lead to very similar results for the S- and P-wave irN 
threshold parameters [111 ] and, consequently, for the coefficients b and d. A similar observation is also made for the 
subleading 27r-exchange NN potential in Ref. 11]. Another interesting result is that both EFTs with and without A 
lead to very similar values for the coefficients a, b and d at NNLO. This seems to contradict the conclusion of Ref. [n| 
where a significant overestimation of the A-contribution to the coefficients b and d by about 25% was found in EFT 
without explicit A fields. As shown in Table HI the coefficients b and d are indeed overestimated in the A-less theory, 
but only by about 7%. There are several reasons for this difference. First, our results for b and d involve contributions 
beyond the ones generated by the A excitation. Secondly, we also take into account the subleading A-contribution 
governed by the 63 + &8 _v ertex from ^j^A wmcn is n °t considered in Ref. (l9| . Switching off the 63 + 6s-interaction, 
however, tend to further decrease the above mentioned overestimation and can, therefore, not explain the discrepancy. 



2 Due to parity invariance, the subleading nNN vertices contain two derivatives or one insertion of M 2 more than the leading one and 
thus need not be considered at NNLO. 
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The most important difference between our work and the one of Ref. [l9| concerns the determination of the LECs C3 
and C4. In [L9J, this was achieved via matching the P-wave irN threshold parameters for j = 3/2, and a^ + . This 
leads to values for 6 and d which are quoted in Eq. (15) of that work. On the other hand, in [ll| we used not only 
the j = 3/2 but also the j = 1/2 P-wave parameters and as well as the S-wave coefficients and 6q, to 
pin down the LECs Cj and 63 + 68. This turns out to be the main reason for the observed difference. 



III. LEADING ISOSPIN-BREAKING EFFECTS 



Our next goal is to study the leading isospin-breaking contributions to the 3NF. This can be done following the lines 
of Refs. I2Q, l2li where we worked out the leading and subleading isospin-breaking 2N and 3N potentials, see also 
Refs. [22I I23L \24. 125I] for a related work on this subject. Here and in what follows, we adopt the same counting rules 
for the electric charge e and the quark mass ratio e = (m u — md)/{m u + ra c i) ~ — 1/3 as in Refs. (20l. [2l|. namely 

£ ~ 6 ~A' (4^) 2 ~A 4 - [6 j 

The leading isospin-breaking vertices have the dimension Aj = 2 and correspond to the charged-to-neutral pion and 
proton-to-neutron mass differences [201 ] and, as will be shown below, mass splittings between the different charge states 
of the A. It is easy to verify using Eqs. (|2.1[) and (|2.2[) that the dominant isospin-breaking 3NF which involves an 
intermediate delta excitation results at order v = 4 from isospin-breaking pion-, nucleon- and delta-mass insertions in 
the first diagram in Fig. [T] Notice that the leading isospin-breaking ititNN vertex of dimension Aj = 2 also generates 
the 3NF contribution at order v = 4 [HIHH which, however, does not involve A-excitations and is, therefore, irrelevant 
for the present study. Thus, in order to proceed with the calculation of the leading isospin-breaking 3NF contributions 
due to intermediate A-excitations, we first need to analyze the delta mass splittings in chiral effective field theory. 



A. Delta mass splittings in chiral effective field theory 



To analyze the delta mass splitting in chiral effective field theory, we need to include hard virtual photons in the 
effective pion-nucleon-delta Lagrangian. Specifically, we are interested in the leading-order virtual photon effects that 
show up as loca l op erators of dimension Aj = 3. To construct these terms, we employ the standard spurion method, 
see e.g. [2(1 H3, Ha 29]. The resulting Lagrangian reads 



ft S tJ (Ql - Ql) + / 2 A Sa (Q+)Q+ + / 3 A % {Q+f 



+/ 4 A (t 1 Q+)(tJQ + ) + ft (t 1 Q-)(t1Q_) 



rpV 
9p.V 1 j 



(3.2) 



where (. . .) denotes the isospin trace, the Q± are defined as in [27| and is the pion decay constant in the chiral 
limit. The factor F% ensures that the LECs / 2 A have the same dimension as the corresponding strong LECs [27l ]. 
Notice that while all couplings and masses appearing in the effective Lagrangian should, strictly speaking, be taken 
at their chiral limit values, to the accuracy we are working, we can use their pertinent physical values. We further 
emphasize that the expected natural size of the LECs / A is 



F 2 f • 



(4tt) 2 



(3.3) 



so that, according to Eqs. (|2.2[) and (|3.1[) . the corresponding vertices are indeed of the dimension Aj = 3. The 
electromagnetic mass term of the delta is readily deduced from Eq. (|3.2[) by considering the terms without pion fields, 
e.g. 



C 



Ai=3 
A7, mass 



-F 



(A A + / 3 A ) e 2 % + f2 TT^ + r 3 ) Stj + ft e 2 S l3 5 



'J3 



9^v Tj 



(3.4) 
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Note that the first term in this equation leads to an overall mass shift in the delta quartet, while the other two 
contribute to various splittings. The relevant isospin conserving and strong isospin violating terms at leading order 
are given by [1(3] 

= -ft 4 (X+ - <X+» kj 9»» Tjf = -ft 4 2Ml e r 3 % Tjf + . . . (3.5) 

where ellipses refer to terms which contain pion fields. Note that in Ref. [l(| the LEC was denoted as 05. Combining 
Eqs. (13. 4p and (13. 5|) . we arrive at the leading strong and electromagnetic isospin-breaking contributions to the delta 
mass term 



£ A, mass, IV— T-f 

with 



1 3 

-8m\ - t 3 6ij - Sm A ~5 i3 S j3 



9nv Tjf , (3.6) 



5m\ = -AM%e<£-F*<?f£ , 

Sml = -^ 2 e 2 / 4 A - (3-7) 

Notice that according to the counting rules in Eq. (|3.ip . the dominant term in Eq. (|3.6p arises from strong isospin 
breaking while the leading electromagnetic contribution is shifted one order higher. At least in the nucleon sector, 
this pattern is in a reasonable agreement with the data and the strong shift gives indeed numerically the dominant 
contribution. Since we are interested here only in the leading isospin-breaking effects due to explicit A degrees of 
freedom, it is, strictly speaking, sufficient to keep only the strong delta mass shift. We will, however, keep both the 
leading strong and electromagnetic delta mass shifts in what follows and thus generate a portion of the subleading 
isospin-breaking contributions to the 3NF. At the order we are working, isospin-breaking contribution to the delta 
self-energy results from the tree diagram with a single insertion from Eq. (|3.6|) . The delta mass splittings can, 
therefore, be directly read off from Eq. (|3.6[) . Switching from the Rarita-Schwinger representation to physical delta 
fields (A++, A + , A , A~) one obtains: 

1 r 1 

m A ++ = m A + - om A , 

1 r 1 1 c 2 

m A + = m A + - dm A + -0m A , 

2 

1 r 1 1 c 2 

m A o = toa - - om A + -om A , 
2 

m A - = toa — — Sm A . (3-8) 

Notice that the isospin- invariant delta mass shift Sm A , defined according to toa =fiA +<5toa, can, in principle, 
be extracted using the chiral limit of the NA splitting Ao ~ 330 MeV from Ref. [3(| and the value of the nucleon 

mass in the chiral limit tojv— 890 MeV from Ref. [3l|. This leads to toa— 1220 MeV. We further emphasize that 
in the absence of electromagnetic corrections (no 5m A -term) , there is an equal spacing between the members of the 
quartet. In that case, one has fh A — rh A and all splittings are equal to Sm A /3. We also recover the SU(6) results 
toa++ — rn A - = 3(toa+ — rn A o) independent of the strength of the electromagnetic LEC f A , see e.g. [3^ . Our results 
for the strong splittings are, of course, in agreement with the ones of Ref. [53]. In that paper, higher order corrections 
were also evaluated with particular emphasis on the quark mass dependence of the delta splittings to be used as chiral 
extrapolation functions in lattice gauge theory. 

To further analyze the delta mass splittings, we need input values for some of these masses. Astonishingly, the 
available information of these masses is fairly scarce and uncertain, see e.g. the most recent listings of the particle 
data group [33| . From pion-nucleon scattering, one can extract toa++ and toao as well as the average delta mass fh A 

toa = - {m A ++ + m A + + m A o + m A -) — m A + -Sm A . (3.9) 

The values for the Breit-Wigner masses toa++ and toa<> quoted in Ref. [33| and based on the determinations from 
Refs. [13, [H, [H, H3] are all in a reasonable agreement with each other (within the given error bars). The only 
exception is the analysis of Ref. [38| . which yields significantly different values for toa++ and toao (but has also 
the largest error bars). From photoproduction reactions, one can, in principle, determine the mass of A + . The two 
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available results for A + quoted in [H| differ, however, significantly from each other. We, therefore, refrain from using 
m A + in our study. Finally, no experimental information is available for m A -. To pin down the values for rh A , 5m a 
and Sm A we proceed in two different ways. First, we use as an input the available data on m A ++, m A o (except the 
values from Ref. [38|) and the average delta mass rh A . For the latter, we adopt the value fh A = 1233 MeV which is 
consistent with the estimation of Ref. [33|, fh A = 1231 . . . 1233 MeV as well as with the most recent determination 
from Ref. 39], m A = 1233.4 ± 0.4 MeV. This leads to 

m A = 1233.4 ±0.7 MeV, 5m A — —5.3 ± 2.0 MeV , 5m A = —1.7 ± 2.7 MeV , (3.10) 

where the error bars result from using different input values for m A ++ and m A o, see Table HU Notice that the 



771 A 


5m A 


8m\ 


m A+ + 


m A+ 


m A o 




input 


1233.63 


-6.15 


-2.50 


1230.55* 


1231.35 


1233.40* 


1236.70 


[34] 


1234.10 


-7.20 


-4.40 


1230.50* 


1230.70 


1233.10* 


1237.70 


[35] 


1233.15 


-4.50 


-0.60 


1230.90* 


1232.10 


1233.60* 


1235.40 


[36] 


1232.75 


-3.30 


1.00 


1231.10* 


1232.70 


1233.80* 


1234.40 


[3JJ 



TABLE II: Delta masses and LECs for various input values of m A ++ and m A o as indicated by the star. Additional input is 
the average delta mass m A = 1233 MeV. All values are given in units of MeV. 

obtained results for Sm A and 5m 2 A are of a natural size. Indeed, based on the naive dimensional analysis, one expects 
\5m\\ ~ \eM%/M p \ ~ 8 MeV and \Sm 2 A \ ~ e 2 M p /(47r) 2 ~ 1.5 MeV. 

As an alternative, one can use in the determination of m A , Sm A and 5m\ the quark model relation [32| 

m A + — m A o = m p — m n (3-11) 

instead of the average delta mass m A . This fixes the value of 5m A and leads to less uncertain results for rh A and 
Sm A , see Table Uni for more details: 



m A 


8m\ 


m A++ 


m A+ 


m A o 


m A - 


ra A 


input 


1232.49 


0.53 


1230.55* 


1232.11 


1233.40* 


1234.43 


1232.62 


[3J 


1232.44 


0.03 


1230.50* 


1231.81 


1233.10* 


1234.38 


1232.45 


[35] 


1232.84 


0.23 


1230.90* 


1232.31 


1233.60* 


1234.78 


1232.90 


[36] 


1233.04 


0.23 


1231.10* 


1232.51 


1233.80* 


1234.98 


1233.10 


[3JJ 



TABLE III: Delta masses and LECs for various input values of m A ++ and m A o as indicated by the star. Additional input is 
the quark model relation (|3.11[) . All values are given in units of MeV. 

m A — 1232.7 ± 0.3 MeV , Sm A — —3.9 MeV , Sm A — 0.3 ± 0.3 MeV . (3.12) 

It is comforting to see that the results of both determinations are compatible with each other. In particular, the value 
obtained for the average delta mass, fh A — 1232.45 . . . 1233.10 MeV, agrees with both the estimation of Ref. [33| and 
the determination of Ref. [3!| . We further emphasize that the LECs and f£ can be deduced from Snrff and (5m™ 1 
if one uses the relation (|3.1ip separately for the strong and electromagnetic mass shifts. Using <5m^ r = (m p — m n ) str = 
-2.05 ± 0.3 MeV and Sm^ = (m p - m„) cm = 0.76 ± 0.3 MeV from Ref. [H, see also [U for a recent determination 
from lattice QCD, one obtains = 3c 5 = -0.24 ± 0.04 GeV^ 1 and / 2 A = 3/ 2 = -2.9 ± 1.1 GeV" 1 , where c 5 and f 2 
are the corresponding LECs in the nucleon sector. 



B. Field redefinitions and the proton-to-neutron mass shift 



Having worked out the delta mass splitting in chiral EFT, it is, in principle, a straightforward task to calculate 
the dominant isospin-breaking 3NF contribution due to explicit A following the line of Ref. [2p| . The calculations 
may be facilitated if one eliminates the nucleon mass shift from the effective Lagrangian. This allows to get rid of 
diagrams which involve reducible topologies and, therefore, enables to use the Feynman graph technique to derive 
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the 3NF. As already demonstrated in Ref. [25| for EFT without explicit A, the proton-to-neutron mass difference 
can be eliminated from £ c g via a suitable redefinition of the pion and nucleon fields (or, equivalently, via a suitable 
local SU(2)y transformation) in favor of new vertices proportional to Smpf. This approach can be straightforwardly 
generalized to include the A as an explicit degree of freedom. To be specific, consider the effective chiral Lagrangian 
£ e ff($, J), where $ = {U, N, N, T, T} (J = {r M , l^, s,p}) collectively denote the pion, nucleon and delta fields (right- 
and left-handed, scalar and pseudoscalar external sources). The effective Lagrangian £ c ff(<£>, J) is invariant under 
local chiral rotations G = SU{2) L x SU(2) R : 



£ eff ($, J) £ cff ($', J') = £ cS ( 5 ($), g(J)) = £ eff ($, J) . 



(3.13) 



For the purpose of computing S-matrix elements in the few-nucleon sector, the external sources can be set to zero 
from the beginning (with the exception of the scalar source that is equal to the quark mass matrix in this limit). In 
the absence of external sources, the effective Lagrangian £ c ff($) = £ c ff($, J)| J=Q is only invariant under global chiral 
rotations. Local chiral transformations g € G will, in general, affect its form 

C cB (*) ^ Arff (<?($)) - £eff($) + 5£ eS {$) , (3.14) 

and obviously may be viewed as a redefinition of fields The task is now to choose the transformation V in such a 
way that the resulting correction 5£ c jf (<f>) eliminates the nucleon mass shift term 

1 



N6m N ^T 3 N, 



5m N = -4c 5 eMZ-e'Ftf 2 



(3.15) 



in £ e ff(<l>). Notice that in practice, it is more convenient to compute the correction 5£ c s applying the inverse rotation 
to the external currents 



£ eff (<?($)) = £ eff ($, g-\j)) = £ cff ($) + 5£ cS ($) 
Following Ref. [Hj], we choose the local SU{2)y transformation such that 

r 'u = Vr^ + iVd M V\ I' = VI ^ + iVd M V*, s' = VsV\ 



V = VpV ] 



(3.16) 



(3.17) 



with V = exp(i v ■ x Sm^r 3 /2). Notice that as pointed out in Ref. (2f|, this transformation does not lead to explicitly 
x-dependent vertices since the interactions in £ e ff conserve electric charge. One obtains for the chiral vielbcin = 



— ir^)u — u(d IJl — iZ^)tt'l and connection T 



1/2 [u^{d^ - ir^u + u(d^ - ilfjvt 



j=o 



j=o 



,7=0 



J=0 



J=0 



J=0 



5m n 
2 

drriN 



(llV 3 H — UT 3 v) 



V^T X 7T] 3 + C(7T 3 ) , 



= r, 



.Stun 



— i- 
j=o 4 

. 5m 

— i 

,/=o 2 



1> M (llV 3 ll + HT 3 M^ 



at o . 5m^ 2 o 

-*V T + ^"CF T J - 7T • T7T 3 ] 



(3.18) 

0(7T 4 ), 



and, consequently, for the nucleon and delta kinetic energy terms: 

Niv-DN = N (iv ■ d + iv -T) N —> N (iv ■ d + _I^_ r 3 + q(tt 2 ) ) X . 
(iv Dij) 9llv TV = -T?(i{vd + vY)5 ij +e ijk {T k vY))g llv T% 



_rpH 



■d 



35m 



N 3 



5 ij+ 0(^) \g^Tjf. 



Here we have used the following identity for the Rarita-Schwinger field: 



if?5 t] T k T» 



Thus, we finally end up with the modification of the lowest-order effective Lagrangian 



■.Ai=0 
-off 



-Ai=0 



5m N [(v ■ dn) x 7r] 3 + ^(5m N ) 2 (n' 2 - 7r 3 7r 3 ) + C(7r 4 ) + 
^Lft(r 3 5 ij + 0(n 2 )) g^Tf. 



SmN 



N (r 3 + 0(tt 2 )) N 



(3.19) 
(3.20) 



(3.21) 
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For the pion and pion-nucleon sectors, these results agree with the ones obtained in Ref. 2£ 



Finally, we would like to emphasize that although the applied field redefinition is, strictly speaking, of a more general 
type than the one discussed in [42|, , it is straightforward to show that both £ e g and £ e ff + SC e g lead to the same 
S'-matrix elements in the few-nucleon sector. The corresponding Green functions in the original and modified theories 
are, however, not the same 3 and related to each other by the transformations V: 



5 2n Z[J, V ,fj] 



<%1 fal) • ■ ■ Sr )jn {Xn)5r)ii (yi) ■ ■ ■ (Vn) 



E 



5 2n Z[J',r l ' ) fj'] 



x V kujl {x x ) . . . V knJn (x n )vl h ( yi )... VI K (y n ) . 

Here Z[ J, 77,77] is the generating functional 

e iZlJ,r,,fj] = [ e t / d l x{C ct{ (®,J)+fiN+Nri\ ^ 



J'—r}'—r}'—0 

(3.22) 



(3.23) 



and 77' = Vi], rf = rjV'. 



C. Isospin-breaking 3NFs due to explicit deltas 



We are now in the position to present our results for the leading isospin-breaking contributions to the 3NF due to 
an intermediate A-excitation. A straightforward evaluation of Feynman diagrams (b-d) in Fig. [T] yields the following 
charge-symmetry conserving (CSC), i.e. class-II in the notation of Ref. [2lj . and and charge-symmetry breaking (CSB), 
i.e. class-Ill, contributions to the 3NF: 



V, 



Class-II 
3N 



E 



9 A h A 



q% ■ cnqk ■ &k 



SMl 



r l8F*A [ql + MX+M* 

2 



-3^3 



ASMl 



A 



q% ■ q k 



V, 



Class-Ill 
3N 



% 

E 

wine 



3Sm 'A W x .,33 

8A [ 1 kl j 



g 2 A h 2 A (5m\-35m N ) q} ■ Siq k ■ a k 



216 J FiA 2 



[q^+MX+M^] 



f* 4A 
<fi x q k ] ■ (jj 

5 [t, x T k ] 3 [qi x q k ] ■ aj 



4(r,- ■ T k T? - 2t { ■ r fc T 3 ) <fi • q k 



(3.24) 



The above expressions together with Eqs. (46), (49), (52) and (54) of Ref. (H) corresponding to the contributions 
from 27r-exchange diagrams without intermediate delta-excitation provide the leading isospin-breaking 3NF in EFT 
with explicit A. 



The obtained results can be straightforwardly transformed into configuration space: 



Class-II 



3N 



E 



9A n A 



2 hi Ml 



2887r 2 FiA 



2- (Ti ■ Vy <7fc ' V k j 



(25MI 
\ Ml 



Ui{x i0 ) U 2 {x kj ) 



4A 



Uifaj) Ui(x kj ) 



(5Ml 
\2Mi 



'Ti X T,] 3 T? 



k Ui(xij) U 2 (x kj ) 



35m 



8A 



[Ti X T k ] 3 T 3 Uxfaj) Ul(x k j) 



Class-Ill 
3N 



E 



g\h 2 A {5m\ - 35m N )M% _ 



34567r 2 FiA 2 



Oi ■ V,j <7fc ■ Vfej 5[t; X Tk] [Vjj X Vjy] • CTj 



4(Tj • T k r 3 - 2r 4 • T fc T 3 ) Vjj • Vfcj Ui(xij) Ux(x kj ) , 



(3.25) 



For example, the free nucleon propagator in the modified theory does not contain 5m jv any more. 
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= -o 



FIG. 2: Diagrams in the A-less EFT which reproduce the A-contributions to the leading isospin-breaking 3NF via resonance 
saturation of the LECs. Filled circles denote isospin-breaking pion, nucleon and delta mass shifts. Filled square and diamond 
refer to isospin-conserving vertices of dimension A, = 1 and Ai = 2. Crossed diamond denotes isospin-breaking strong and 
electromagnetic vertices of dimensions A, = 4 and A, = 5, respectively. For diagrams involving the insertion of the nucleon 
mass shift, only two representative graphs are shown. 



where = f. \ — r) is the distance between the nucleons i and j, x*i = M^ri, Vi act on xt and Xij = Further, 
we have introduced the profile functions 



Air f d 3 q e^ s / M ^ 



e 



(2tt) 3 q 2 + Ml x 



w = ii ^ l(^WTMW = e ~ X - (3 - 26) 



It is instructive to understand how the obtained 3NF contributions are reproduced in EFT without explicit deltas. The 
trivial A -1 - and A _2 -dependence of the 3NF on the delta-nucleon mass splitting arising from the static propagator of 
the delta which enters Feynman diagrams (b-d) in Fig. [T] ensures that the results are reproduced in the A-less theory 
by a finite number of higher-order graphs (via resonance saturation of certain LECs). This is depicted in Fig.[2l Here, 
we have switched back to the original Lagrangian with the nucleon mass shift which is more convenient to discuss 
resonance saturation. It should, however, be understood that for Feynman diagrams involving the nucleon mass shift, 
only irreducible contributions are taken into account. Consider first the CSC 3NF proportional to SM 2 in Eq. (|3.24|) . 
As expected, this contribution is exactly reproduced in the A-less theory by an appropriate shift in the LECs C3 and 
C4 [T4I, C3 = — 2c4 = —Ah A /(9A), which in this case is the subleading 4 effect (i.e. order v — 5). The remaining 
terms in Eq. (|3.24[) are proportional to A~ 2 which means that they arise in the A-less theory only at sub-subleading 
order v = 6. This is consistent with the absence of CSB contributions proportional to 03,4 Sttin at order v = 5, see 
[II], [44)]. It is easy to verify that terms in Eq. (|3.24[1 proportional to Stun are indeed reproduced in the A-less theory 
by A-saturation of the isospin-conserving sub-subleading pion-nucleon vertices (a complete list of these terms in the 
Lagrangian can be found in Ref. [45| ) while the ones proportional to 5m\ and 8m\ arise from resonance saturation 
of the sub-subleading strong and electromagnetic isospin-breaking pion-nucleon vertices, see FigfJl 

To get a rough idea about the size of the isospin breaking 3NF contributions to e.g. the 3N binding energy, one 
can look at the strength of the corresponding r-space potentials in Eq. (|3.25|) . For the CSC terms oc 5M 2 one gets 
5M 2 g 2 A h 2 A M^j '(1447r 2 F^A) ~ 50 keV. Numerically, this is expected to be the biggest isospin-breaking 3NF effect. 
Notice that in the theory without explicit deltas, this contribution is shifted to the subleading order v = 5 [21[ ■ The 
strength of the remaining CSC 3NF contribution cx 8m\ is much smaller, |<5m^| g 2 1 /i 2 1 M^/(3847r 2 F 4 A 2 ) ~ 1.5 keV 
(here we use the central value for 5m 2 A = -1.7 MeV from Eq. ((3~Tf])) V The estimated size of the CSB 3NF in Eq. ([3T23)) 



4 Numerically, however, it is expected to give the strongest isospin-breaking 3NF [21|. see also the discussion in the next paragraph. 
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using 8m\ = —5.3 MeV is 

\Srni-36m N \ ^J F £ 2 ~3keV, (3.27) 

which is comparable to the typical size of the leading CSB 3NF obtained in Ref. [2l[ and based on EFT without 
explicit deltas, 5m N g\M*/(256ir 2 F£) ~ 7 keV. We further emphasize that the CSB 3NF in Eq. ([3~2"4"j) vanishes 
exactly if one adopts the quark model relation (|3 . . 
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